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Abstract 

Much of the nontrivial dynamics of the one dimensional Complex Ginzburg-Landau 
Equation (CGLE) is dominated by propagating structures that are characterized 
by local "twists" of the phase-field. I give a brief overview of the most important 
properties of these various structures, formulate a number of experimental challenges 
and address the question how such structures may be identified in experimental 
space-time data sets. 
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1 Introduction 



A large body of theoretical work lias been devoted to unravel the intricate 
space-time dynamics of the ID CGLE. One of the most exciting developments 
over the past few years has been the understanding of various aspects of the 
non-trivial, fully nonlinear behavior of this model in terms of the properties 
of a number of "coherent structures". While much earlier work has focused 
on so-called Nozaki-Bekki (NB) holes [1], two other, recently revealed fami- 
lies of coherent structures appear to play the dominant role in large parts of 
parameter space. These structures are the weakly nonlinear Modulated Am- 
plitude Waves (MAWs) that occur when plane waves become linearly instable 
[2,3] and the related, but more nonlinear, Homoclons which are connected 
to defectQ formation [4-6] . The work on these two structures has been pub- 
lished rather recently and it is therefore no surprise that, at present, their 
experimental relevance is unclear. 



^ The defects that occur in the ID CGLE are sometimes referred to as phase-slips; 
these defects do not persist, but rather occur at isolated points in space-time. 
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Coherent structures have a fixed spatial profile, and their dynamics is a combi- 
nation of propagation and oscillation; in CGLE-language, A{x, t) = exp*~"'=°''* A{x— 
Vcoht)- One would expect that coherent structures observable in experiments 
are (i) linearly stable (such that they are attracting) and (ii) structurally sta- 
ble (such that small perturbations of the equations of motion do not destroy 
them). If a certain coherent structure is (structurally) unstable, one expects to 
see, in experiments or simulations, that the "shape" of the structure changes 
over time: I will refer to such structures as incoherent structures. One may 
require that by a smooth change of parameters or initial conditions, a certain 
incoherent structure can be brought arbitrarily close to its coherent counter- 
partQ; the slower the nontrivial dynamics is, the "more coherent" a certain 
structure is. 

The crucial complication one encounters when confronting theoretical predi- 
cations for the behavior of local structures to real experimental data is that 
Homoclons and MAWs are (almost) always linearly unstable [2,3] (see also 
section 2.2.1), while almost all NB holes are structurally unstable [7]. Never- 
theless, the unstable structures are important building blocks for the dynamics 
of the CGLE [2-6]; in many states incoherent structures occur, and these are 
often related to the unstable coherent Homoclons and MAWs. Therefore the 
study of experimental space-time data sets is, I believe, essential for the un- 
derstanding of ID wave systems: snapshots of the field simply do not contain 
enough information [8]. In that sense, the situation in one dimension is more 
difficult than in two dimensions, where a central role is played by spirals that 
can easily identified in snapshots of the field [9]. 

Some additional problems one encounters when comparing experimental data 
to theory are: (i) Many of the theoretical studies have focussed on the chaotic 
regimes, giving the false impression that MAWs and Homoclons only play a 
role in chaotic states, (ii) In most experiments the values of the linear and 
nonlinear dispersion coefficients, which are essential in the theoretical descrip- 
tion, are not known. (Hi) There is some confusion, I believe, concerning the 
relevance of so-called Nozaki-Bekki holes [1] . Since their analytical form has 
been known for more than 15 years, these holes have been studied extensively 
[7,10,11] and there are some claims in the literature that these are observed 
in experimental systems [12-14]; as I will argue below, it may be beneficial 
to have a second look at some of the data. To clarify this situation, I will 
suggest a number of experiments designed to probe the relevance of MAWs 
and Homoclons in traveling wave systems. In addition I will discuss how to 
distinguish incoherent Homoclons and NB holes. 

The outline of the paper is as follows. In section 2 a brief theoretical intro- 



^ this is more difficult for structural instabilities, which are due to perturbations of 
the underlying dynamical system that may not so easily be reversed 
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duction to the CGLE and the coherent structures framework is given, and the 
properties of the following coherent structures and ODE orbits are discussed: 
(i) Plane waves (corresponding to fixed points) (ii) Homoclons (correspond- 
ing to homochnic orbits) (iii) Modulated Amplitude Waves (MAWs) (limit 
cycles) (iv) Nozaki-Bekki (NB) Holes (heteroclinic orbits) . In Section 3 I give 
an overview of the dynamical behavior of the incoherent MAWs and Homo- 
clons and suggest a number of experiments to probe their properties. Section 
4 contains conclusions and a short outlook. 



2 The ID CGLE 

2.1 Basic properties 

The one- dimensional complex Ginzburg Landau equation describes pattern 
formation near a supercritical Hopf bifurcation [11,15]. The focus of this pa- 
per is on experiments that produce a single ID traveling wave (or a uniform 
oscillation) via a forward Hopf bifurcation. In many of such wave systems, 
states with both left and right traveling waves occur, and an overview of the 
behavior of the coupled Ginzburg-Landau equations that describe this system 
can be found in [16]. Here it is supposed that the system can be manipulated 
so as to contain a single travehng wave state, such that the ID CGLE is the 
appropriate amplitude equation. In its full dimensional form this equation 
contains a large number of coefficients, most of which can be scaled out for a 
theoretical analysis (see appendix A). After such a rescaling, the CGLE reads: 

dtA^A+{l + ici)d^^A-{l-zc3)\A\^A . (1) 

The CGLE displays a wide range of behavior as function of the coefficients 
Ci and C3. For example, when ci and C3 are of opposite sign, the dynamics 
of the CGLE is essentially relaxational, while when ci and C3 go to 00, it is 
integrable [11,15]. Away from these limits, the dynamics interpolates between 
ordered and chaotic [18,19]. 

An important symmetry of the CGLE is its phase invariance; if A is a solution 
to the CGLE, so is Ae"^ (for constant 0). This symmetry is related to invari- 
ance of the underlying system with respect to shifts in (space-) time. Writing 
A in its polar representation as A{x,t) = a{x,t)e^'^^^''^\ only the derivatives of 
the complex phase are relevant, and it is helpful to think in terms of the mod- 
ulus a and the phase-gradient dx4>. This latter phase-gradient is also referred 
to as a "local wavenumber" , denoted by q. 
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The simplest nontrivial solutions to the CGLE are plane waves of the form 



A= ^l-ql^ exp{i{qpyjX - uJpyjt)), uJpw = -C3 + ^^^(ci + C3) (2) 

Note that the local wavemimber of plane waves (2) is constant: q{x,t) = 
These states are the background for the dynamics that will described here. 
A hnear stabihty analysis reveals that these waves are prone to the Eckhaus 
instability when [20-22] 



^p- > 3 _ C1C3 + 24 • 



The band of wavenumbers for which plane waves are stable is widest for ci = 
C3 — O and shrinks when these coefficients are increased; when C1C3 > 1 there 
are no linearly stable waves left and chaos occurs. 

The importance of MAWs and Homoclons can be understood by consider- 
ing the dynamical fate of local perturbations of the background wavenum- 
ber of a plain wave consisting of a "twist" of the A field of the CGLE, or, 
equivalently, a local concentration of q. Inside the stable band, the linear evo- 
lution of the local wavenumber is given by a combination of diffusion and 
advection [11]: qt = Dq^x + Vgrqx, where Vgr is the non/mear group- velocity: 
Vgr = dupw/dqpw = 2(ci + 03)5^^. Such linear analysis cannot capture the case 
of nonlinear phase-twists, nor the case when the phase-diffusion coefficient D 
becomes negative (which happens outside the stable band). As will be dis- 
cussed below in more detail, the general evolution of phase-twists is, to a large 
extend, governed by the existence of the MAW and Homoclon coherent struc- 
tures as illustrated in Fig. 1 [2-6,21-23]. (i) For a wide range of parameters, 
both when the background wave is stable or unstable, there exists a nonlinear 
coherent structure called a "Homoclon' that corresponds to a particular local 
phase-twist structure. This structure is linearly unstable, and acts as a scpa- 
ratrix: "smaller" phase-windings decay, while "larger" ones evolve to defects 
(arrows A and B in Fig. 1). (ii) There exist closely related but less nonlin- 
ear structures referred to as Modulated Amplitude Waves (MAWs) . For small 
background wavenumber, the MAWs occur via the forward Hopf bifurcation 
that occurs when a plane wave turns unstable; in this case, the instability 
of the laminar background carries over to the MAWs, and MAWs are often 
linearly unstable [2]. Unstable MAWs often lead to a disordered state called 
phase-chaos, in which patches of transient MAWs occur [2]; the structure of the 
phase-gradient peaks in phase-chaos is comparable to those of MAWs (arrows 
B and C in Fig. 1). For large background wavenumber, MAWs can occur via a 
subcritical bifurcation, and may become linear stable [2,3,23]. (in) Homoclons 
and MAWs are closely related and disappear in a saddle node bifurcation for 
sufficiently large Ci and C3. After this has happened, arbitrarily small phase- 
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Fig. 1. Schematic diagram illustrating the relation between MAWs, Homoclons and 
CGLE dynamics for a background wave of small wavenumber. The vertical axis can 
be thought of as the maximum phase gradient of a local structure, while the hori- 
zontal axis represents the coefHcients ci and C3. The curves represent stable plane 
waves (continuous line) , unstable plane waves (dashed line) , "Upper branch" struc- 
tures (dot-dashed curve) and "lower branch" structures (dotted curve), while the 
arrows represent typical incoherent dynamics (see text). Homoclinic and limit-cycle 
structures (i.e., isolated holes and periodic modulated structures) occur on both 
branches, but the most relevant of these are MAWs, periodic structures on the lower 
branch, and Homoclons, isolated structures on the upper branch. Isolated structures 
on the lower branch can be called P ^ 00 MAWs, while periodic structures on the 
upper branch can be called "periodic Homoclons"; neither play an important role 
in dynamical states. For more details see [2-6] 

gradients diverge (arrow D in Fig. 1) and defects form spontaneously [2-5]. 

MAWs and Homoclons both are characterized by local concentrations of phase 
gradient and corresponding dips of \A\. To complicate matters, there is an- 
other coherent structure, known by the name Nozaki-Bekki hole, which is also 
characterized by a dip of A, and which is a source of waves with different 
wavenumber. Before turning our attention to the experimental relevance of 
these structures, a brief overview of the coherent structures framework and 
the properties of these three families of coherent structures is given. Readers 
familiar with these structures can skip this introduction and go straight to 
section 3. 
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2.2 Coherent structures 



In this section I will briefly discuss the coherent structure framework and list 
the most important properties of coherent MAWs, Homoclons and NB Holes. 
Section 3 focuses then on the incoherent structures and their relevance for 
experiments. 

The temporal evolution of coherent structures in the CGLE amounts to a 
uniform propagation with velocity Vcoh and an overall phase oscillation with 
frequency uJcoh [24] : 

A{x, t) = e-'^-^'a{i)e''^^^\ ^ :^ x - W (4) 



When the ansatz (4) is substituted into the CGLE, one obtains a set of 3 
coupled first order real ordinary differential equations (ODE's) (see appendix 
B). These equations can be written in a number of forms; the representation 
used here employs a, the local wavenumber q := d^(f), and k := {l/a)d^a as 
dependent variables. Orbits of the ODE's correspond to coherent structures 
of the CGLE. 

The ODE's (B.1,B.2) allow for a number of fixed points. For fixed Vcoh and 
u!coh, there are two fixed points with a ^ 0, and these correspond to plane 
waves: 



Qpw = {^yv'^oh + 4:{^coh + C3)(Ci + Cs) + Vcofe)/2(ci + C3) (5) 

a=^^-Qpw (6) 
K = (7) 

The relation between these two fixed points is that their corresponding plane 
waves have the same frequency, Ucoh, in the frame moving with velocity Vcoh- To 
see this, note that plane waves in the coherent structures framework (4) are of 
the form: exp{—iu!coht)0'{0 exp(i0(x — feo^jt)), where = QpwX^- Comparing 
this to a plane wave oc exp^igp^x — Up^t) one finds that in the stationary frame 
the frequencies of the plane waves given by Eq. (5)-(7) are: Upw — cUcoh+QpwVcoh- 
Demanding that the dispersion relation for plane waves given by Eq. (2) is 
satisfied yields: Up^, = lUcoh + VcohQpw = -C3 + gp^(ci + C3), from which Eq. (5) 
immediately follows. The role of the parameter ujcoh in the coherent structures 
ansatz can therefore be interpreted as follows. If cUcoh would be 0, than the 
phase velocity of the plane waves would equal the propagation velocity of 
the coherent structure. Due to the phase-symmetry of the CGLE, these two 
velocities may (and usually will) be different, and in such case cucoh 7^ 0. 

Below, three famihes of coherent structures are discussed: MAWs (correspond- 
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ing to limit-cycles), Homoclons (corresponding to homoclinic orbits) and NB 
holes (corresponding to heteroclinic orbits). 



2.2.1 Limit- cycles and MAW s 

Limit cycles are periodic solutions to the ODE's (B.1,B.2). These orbits are 
structurally stable and generically persist when the ODE's are perturbed. In 
particular, once a limit cycle is obtained for certain values of uJcoh and Vcoh-, 
limit cycles generically exist for nearby values of ijJcoh and Vcoh- For many 
parameter regions, there are two distinct limit-cycles that can be identified 
by their "size" in phase space; for fixed oJcoh-, small(large) orbits occur for 
small(large) values of v^oh (see Fig. 2). The smallest of these orbits occur when 
one of the plane wave fixed points of the ODE's undergoes a Hopf bifurcation 
(and possibly a subsequent drift-pitchfork bifurcations, see [2,3]). A weakly 
nonlinear analysis shows that this bifurcation is supercritical (forward) only 
when cf (1 — 6c|) + Ci(2c| + 16c3) — (8 + C3) > 0; in many cases, this bifurcation is 
subcritical [21,3]. In the CGLE, these small cycles correspond to periodically 
modulated plane waves that we refer to as MAWs. They can be characterized 
by their spatial period P and the "average winding-number" i/, which can be 
defined as 1/P /^f dxdxcj) [2] . When following a particular branch of solutions, 
P and p are functions of Ci, C3, Vcoh and Ucoh and can be obtained by a numerical 
analysis of the ODE's [2,3]. (Fig. 2a). 

The larger limit-cycles correspond to strongly nonlinear modulations of plane 
waves in the CGLE (Fig. 2b). While there is a whole family of these, again 
characterized by P and the limit where the period P goes to infinity is 
most relevant!^ this particular family of CGLE solutions are referred to as 
Homoclons (see below). The two families of limit cycles merge and disappear 
in a saddle-node bifurcation; for details see [2,3]. 

Since there is no closed expression for the full MAW family available, a detailed 
analysis of their range of existence and stability necessarily involves numerical 
calculations. Even when one limits oneself to the "small, weakly nonlinear" 
branch, one still would like to obtain the existence and spectrum as a function 
of Ci, C3, P and z/. Here I will summarize the main results; more details can be 
found in [2,3]. 

• For fixed P and z/, there is a certain range in ci, C3 space where the MAWs do 
exist. For large Ci and C3 this range is limited by the aforementioned saddle- 
node bifurcation, while for small Ci and C3 it is essentially the Hopf/pitchfork 
bifurcation that limits their existence. 



^ Since Homoclons are unstable, finite P solutions rapidly evolve to left and right 
traveling holes, which are well-approximated by interacting, infinite P solutions. 
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Fig. 2. (a) For ci = C3 = 1.1, the ODE flow in a,q space of a "little" limit cycle 
for Vcoh = 0.381 and iOcoh = — l-ll- (b,c) Corresponding |^| and q profiles of the 
coherent MAW characterized by a period P ~ 15 and winding-number v ~ 0.1 (d) 
For ci = C3 = 1.1 one can also obtain a "large" limit cycle for Vcoh = 0.8864 and 
(^coh = — l-l- (e,f) Corresponding \A\ and q profiles of a "periodic Homoclon" with 
P^20,u^ 0.13. 

• For small all MAWs are linearly unstable. There are roughly two in- 
stability mechanisms that compete. For small P, an "interaction" mode is 
dominant. This mode acts as to break the periodicity of the MAWs, but 
the number of phase twists stays the same [2]. For large P, the "split- 
ting" instability is dominant. This latter basically is the Eckhaus instabihty 
acting on the long patches of plane wave that occur for large P; this in- 
stability tends to generate new phase inhomogeneities between the peaks 
of the MAWs, roughly doubling the number of phasetwists. The dynamical 
state that occurs for small u due to the competition of these instabilities is 
called phase chaos, which, for ci and C3 close to the so-called Benjamin-Feir- 
Newell curve (ciCa = 1) can be approximated by the Kuramoto-Sivashinsky 
equation [2,11,15]. I am not aware of any experimental observation of phase 
chaos in one space dimension. 

• For larger values of MAWs may become stable and phase chaos can 
be suppressed. Such stable structures have also been called "compression 
waves" or "wound up phase-chaos" [2,3,21,23,25]. 
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2.2.2 Homoclinic orbits and Homoclons 

Homoclinic orbits start from one of the plane wave fixed points, make an 
excursion througli pliase space and then return to the same fixed point. They 
are structurally stable and occur as a one-parameter family. Once such an orbit 
is obtained for certain values of Vcoh and ujcoh, when Ucoh is varies, Vcoh has to 
be varied accordingly to obtain a homoclinic orbit again. These homoclinic 
orbits can be obtained by letting the period P of the limit cycles diverge. 
Similar to these limit cycles, there are two branches of homoclinic solutions; 
the small ones will be referred to as P — oo MAWs, while the "large, strongly 
nonlinear" of these two correspond to HomoclonsQ. These are propagating 
structures that connect two plane waves of equal wavenumber in the CGLE. 
Their core region is characterized by a concentration of local wavenumber q 
(or equivalently, a twist of the phase 0), and a dip of the amplitude \A\ [4,5]. 
They can be seen as local phase twists that glide through a background plane 
wave; often their total phase twist is of the order of, but not precisely equal 
to, TT. Their winding number u simply corresponds to the wavenumber of the 
asymptotic plane wave. 

For fixed Ci,Cs and z/, there are left and right moving homoclons, related by 
left-right symmetry of the CGLE (which maps q ^ — g). Their propagation 
velocity is typically much larger than the nonlinear group- velocity of the plane 
wave they glide through; they are neither sources nor sinks. In the comoving 
frame, they experience an incoming wave ahead and leave an outcoming wave 
behind. 

For positive ci and C3, solutions occur for positive (negative) Vcoh when the q- 
profile peak is positive (negative). When ci and C3 are both negative it follows 
(by complex conjugation of the CGLE) that the signs of velocity and g-profile 
peak are opposite. 

As for the MAWs, no analytic expression exists for the Homoclons, and a 
numerical analysis has revealed the following key points: 

• The range of existence of Homoclons is limited by the same saddle-node 
bifurcation as the MAWs. Once a Homoclon has been obtained for certain 
values of ci and C3, this solution persists when the coefficients are decreased 
towards zero; in fact these Homoclon-solutions smoothly connect to the 
analytically known saddle-point solutions of the real (ci = C3 = 0) Ginzburg- 
Landau equation [22]. For large values of Ci and C3, no Homoclons exist, 
and generic initial conditions form defects; the CGLE is then in the defect 
chaotic regime [2]. 

^ The strongly nonlinear but periodic structures do not play an important role and 
hence do not have their own name; they can be referred to as "periodic Homoclons" , 
although this is not very elegant. 
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• The spectrum of the Homoclons consists of two parts, a continuous part 
associated with the asymptotic plane waves, and the discrete part associated 
with a few core modes [4]. For aU cases that I am aware of, there is one 
single unstable core mode; this unstable eigenmode is associated with the 
aforementioned saddle-node bifurcation. 

With hind-sight, the existence of Homoclons and their main properties are 
not entirely surprising. Suppose, for some values of ci and C3, that one can 
construct two "generic" initial conditions, one that does not evolve to defects, 
another that does form a defect. When C1C3 < 1 this is certainly possible: 
the state A — 1 is stable and nearby initial conditions will not form a defect, 
while a state where A = 1 except for an interval where A ~ exp{iQx) will 
form defects when Q is sufficiently large. Now imagine that one interpolates 
between these two initial conditions; let's call the interpolation parameter 
b. In the simplest case there will be a single transition value for b where the 
dynamics changes from non-defect to defect forming. Clearly, at the transition 
point something special must happen, and the simplest scenario here is that 
the time it takes for a defect to form diverges there. This may happen in a 
variety of ways, but the simplest case would be the formation of a saddle-point 
like structure: the Homoclon. Since one needs only one parameter to vary, it 
is reasonable to expect a single unstable eigenmode. 

Note that in the defect chaos regime this scenario breaks down, since only non- 
generic initial conditions (perfect plane waves) will not form defects: indeed, 
no Homoclons exist here. 

This scenario can be pTit on firmer ground for the real Ginzburg-Landau equa- 
tion, for which analytical results are available [22]. The dynamics of this equa- 
tion is governed by a Lyapunov functional that is decreasing over the course 
of time. Of course, one can create saddle points that correspond to unstable 
coherent structures. The simplest example of this is the standing hole solu- 
tion A = tanh(yT/2x), which indeed corresponds to a homoclinic orbit of the 
ODE's (B.1,B.2). It is a numerically easy task to trace a coherent Homoclon 
orbit as a function of Ci and C3, and one finds that in the Ci — C3 — limit, 
a, u — Homoclon smoothly deforms to this saddle-point solution of the real 
Ginzburg-Landau equation; the propagation velocity of the homoclons thus 
goes to zero in the relaxational limit. 

2.2.3 Heteroclinic orbits and Nozaki-Bekki Holes 

The heteroclinic orbits that are relevant here start at one plane-wave fixed 
point and end up on another plane-wave fixed point. They can be obtained 
in a closed analytical form [1,24] and the corresponding CGLE structures 
are the so-called Nozaki-Bekki (NB) holes which connect waves of different 
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wavenumbers qi and q2- 



For the unperturbed CGLE, once such an orbit is found for certain values of 
Vcoh and ujcoh, there will generically be heteroclinic orbits for nearby values of 
Vcoh and cucoh- As pointed out in [24], this is not the generic situation one ex- 
pects (on the basis of counting arguments) for these orbits. Indeed it was shown 
that small perturbations of the CGLE destroy most of the NB-orbits/holes [7], 
and one is left then with a discrete family of solutions: almost all NB holes are 
structurally unstable solutions to the ID CGLE. The dynamical states that 
occur in this situation are discussed in [7] . 

Once qi and q2, the wavenumbers of the adjacent plane waves, are fixed, Vcoh 
and LUcoh can be obtained from a simple phase matching argument, which 
states that in the comoving frame the apparent frequency of the two waves 
needs to be equal to cUcoh- Some manipulation yields that v^b — {qi + Q'2)(ci + 
C3) and ujnb = — Q'iQ'2(ci + C3) — C3. Comparing this propagation velocity to 
the nonlinear group velocities of the adjacent plane waves, 2gi(ci + C3) and 
2?2(ci + C3) respectively, it immediately is clear that the propagation velocity 
of the NB holes is in between these group- velocities. A more thorough analysis 
[24] shows that NB holes are always sources, i.e., in the comoving frame they 
send out waves. 

Since the NB holes are known in closed form, a number of results for the range 
of existence and stability are available. In particular, the range of stability of 
NB holes is limited by the instabilities of their adjacent plane waves on one 
side, and a core instability on the other side. This yields a band of stable NB 
holes, which is mainly located in the regime where ci and C3 have opposite 
signs [11]. 



3 Incoherent dynamics and experiments 

When one is in the lucky position that an experiment shows interesting space- 
time dynamics of local structures, some information on how to identify the 
structures that occur can be found in section 3.4. But if the system by itself 
just produces simple plane waves, which often seems to be the case, some ma- 
nipulations need to be done. The good news is, that for all coefficients of the 
CGLE it is possible to generate transient Homoclons by locally perturbing the 
system, and when the wavenumber of the plane states can be controlled, also 
a number of MAW states can be generated. 

In most cases one would expect to observe incoherent rather than coherent 
MAWs, Homoclons and NB-holes. From the wide range of behavior that oc- 
curs in the CGLE, I will highhght a number of states that hopefully will be 
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accessible in experiments. The examples given below are intended to inspire 
particular experiments, and are not fully worked out recipes. Before going on, 
let me briefly discuss some of the properties of candidate systems. 

• Boundary conditions. Boundary conditions will play an important role. 
Some systems, such as Rayleigh-Benard wall-convection in rotating cells 
[26,27] or sidewall convection in annular containers [25] have periodic bound- 
ary conditions. In this case the linear group-velocity term of the CGLE 
(see appendix A) can be removed by going to the comoving frame. The 
main finite size effects to be expected are then the discretization of possible 
wavenumbers and the fact that a source will necessarily be accompanied 
by a sink. The left and right boundaries of long rectangular systems, such 
as heated wire convection cells [28,29] and sidewall convection [14,30], may 
have a more severe effect on the dynamics. First of all, the group- velocity 
can no longer be ignored here. In addition, in some systems the boundaries 
may act as sources that send out waves of a selected wavenumber [16,17]. 
Of course this can be used in experiments, in particular when the selected 
wave can be made linearly unstable. 

• Control-p£trameters. For all systems, the dimensionless distance to thresh- 
old, e (see appendix A) is experimentally accessible. While e may be scaled 
out of the equation for the amplitude A, it does effect the spatial and tem- 
poral scales. For example, a sudden change in e has the effect of virtually 
"stretching" or "compressing" the spatial scale of A. This could be used 
to manipulate the effective wavenumber and the period P of MAWs (see 
below) . 

When the group-velocity cannot be ignored due to boundary effects, 
changes in e can change the instability of plane waves from convective to 
absolute. The amplitude equations no longer scale uniformly with e in this 
case, and the stability of sources pinned at boundaries may change (for more 
on this in the context of left and right travehng waves see [16,17]). 

In some cases, other experimental parameters are accessible, like the depth 
of the heated wire below the surface [28,29]. Such changes will affect the 
coefficients of the CGLE, although often these are not known, let alone their 
dependence on experimental parameters. 

• Local manipulation. Usually the initial conditions cannot be chosen at 
will, in stark contrast to what is done in numerical studies. Nevertheless, 
some manipulations of dynamical states are usually possible. In convection 
systems, a short burst of local heating can be used to perturb the system. 
Such perturbations can be used to generate fronts (see appendix A) or 
possibly to generate Homoclons (see below) . One may also imagine that by 
periodic local heating one can manipulate wavenumbers or the period of 
MAWs. In systems with a free fluid surface, mechanical manipulations such 
as briefly touching the surface or depositing a few drops of oil [25] may also 
be used as (crude) ways of perturbing the system locally. 

• CGLE Coefficients. For most systems, the coefficients of the CGLE have 
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not been calculated, making a detailed comparison to theory rather cum- 
bersome. By careful experiments (see appendix A) one should be able to 
get good estimates for most of their values [14,26,29,31]. 

3.1 Incoherent MAW dynamics 

A comprehensive overview of MAW properties can be found in [2,3] and I will 
mention the most relevant properties in the list of suggestions for experiments 
given below. 

3.1.1 Stable and unstable MA Ws 

When plane waves turn unstable against the Eckhaus instability, MAWs occur. 
When the wavenumber of the initial plane wave is sufficiently large (i.e., for ci 
and C3 far away from the BFN curve), the bifurcation to MAWs is subcritical 
and such MAWs may be hnearly stable [3,21,23]. For such structures one can 
obtain their velocity, period and winding- number (provided the CGLE coef- 
ficients are known), and confront these with numerical calculations of MAW 
profiles. In the subcritical regime, it may be possible in this regime to excite 
MAWs by sufficiently strong local, periodic forcing of the waves. Such exciting 
of MAWs can of course also be done at the boundary of a system, although 
periodic boundaries are probably best suited for such experiments. 

Starting from stable MAWs, a number of MAW properties may be probed. 

(1) When the period of MAWs is altered (by periodic local heating or sud- 
den changes of e) and increased, two possible types of dynamics may 
occur. When ci and C3 are relatively small, the MAW may become prone 
to the splitting instability [2,3], and new peaks in q may be generated 
between the old ones, effectively decreasing the period P. When Ci and 
C3 are larger, the effect of increasing P may be to cross the saddle-node 
bifurcation, and defects will then occur [2] . Since the occurrence of this 
saddle-node appears to be a crucial ingredient for the understanding of 
the transition from phase to defect chaos, this would be extremely inter- 
esting to see. 

(2) When the period of MAWs is decreased, the interaction instability may 
kick in, which breaks the periodicity of the MAWs [2]. 

(3) When the winding number of MAWs becomes sufficiently small, MAWs 
become unstable (see [3,23]). In principle, the effective wavenumber changes 
when e is varied, but at the same time the effective period P changes also. 
Therefore, a change of e can either lead to MAW instabilities or defect- 
formation, depending on which instability is closest. If it is possible to 
manipulate P by local periodic heating, then changing this period ac- 
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Fig. 3. Space-time plot of the evolution of the Benjamin Feir instability of a plane 
wave with initial wavenumbcr qpy^ = 0.41, for ci = C3 = 0.8 in a system with periodic 
boundaries and size 400; only part of the system is shown. The slow growth of a 
transient MAW can be observed, and since this MAW is "beyond the saddle-node" , 
it does not saturate but develops defects. This dynamics is the one depicted by arrow 
D in Fig. 1, and is very reminiscent of defect formation observed in hydrothermal 
waves (see Fig. 6 of [25]) and rotating Rayleigh-Benard convection [26]. 

cording to the changes made in e may be a way to change the effective 
value of v. 

(4) While it is difficult to say anything general about the MAW profiles, it is 
known that the minimal value of 1^41 and the extremum of |g| both grow 
as a function of the period. In fact, near the SN they decrease {\A\) and 
increase {\q\) sharply, which could be used to estimate if this bifurcation 
is close by. 

3.1.2 Transient MAWs. 

When the plane wave is stable but close to its Eckhaus instabihty, perturba- 
tions of this plane wave will decay rather slowly, and transient MAWs can be 
formed. By measuring the decay time as a function of the wavenumber, one 
can make an estimate for the proximity of the Eckhaus instability (see [25]). 

When a plane wave becomes unstable due to the Eckhaus instabihty, the first 
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stages arc characterized by the development of periodic modulations. This 
period is determined by the most unstable wavenumber of the Eckhaus in- 
stability. When these modulations do not grow out to phase chaos or MAWs, 
but instead grow without bound, defects will eventually occur; a typical ex- 
ample is shown in Fig. 3. The transient MAWs are, as usual, characterized 
by their period P and winding-number v (which is equal to the wavenumber 
of the initial wave). Defect formation is then expected to occur when, for the 
current values of Ci and C3, P and v lie beyond the appropriate saddle- node 
bifurcation [2,3]. When P and v are relatively close to the saddle- node, the 
formation of defects can be rather slow, and transient holes can be observed; 
these are clearly incoherent (see [25,27] for possible experimental realizations 
of this) . 



3.2 Incoherent Homoclon dynamics 

Homoclons are never linearly stable; over the course of time, they either slowly 
decay or grow out to form a phase-slip [4,5]. It will be assumed that their 
dynamics is studied in a background of linearly stable plane waves. In earlier 
studies it was found that the instability of the Homoclons is due to a single 
"core"-mode, and its unstable eigenvalue is the one which changes sign at the 
saddle node bifurcation where MAWs and Homoclons merge and disappear 
[2,3]. 

This single weakly unstable eigenmode is reflected in the dynamics of the 
Homoclons. Many sufficiently localized wavenumber-blobs will be attracted 
to the ID unstable manifold of the Homoclon; subsequently they then evolve 
along this manifold, in either the "decay" or the "phase-slip" direction. One 
can loosely think of the homoclinic holes as unstable equilibria, or separatrices, 
between plane waves and phase-slips [4,5]. 

These properties can be probed as follows: Suppose one has a stable plane 
wave, and locally perturbs it by generating a twist in the phase field of A. 
When this twist is small, the linear stability of the plane waves will govern the 
dynamics and such a perturbation will decay diffusively (downward pointing 
arrow A and B of Fig. 1). However, when it is strong enough, and the initial 
conditions "passes through the Homoclon separatrix" , the twist will grow out 
to form a defect (upward pointing arrows A of Fig. 1). When tuning the initial 
twist so as to be in between the decay and defect scenario, one can obtain a 
Homoclon that exists for an arbitrarily long time (In practice, noise may limit 
this time.). The non-trivial prediction is that, since there is only one unstable 
eigenmode, one only needs a one-parameter family of initial conditions to "hit" 
the unstable Homoclon solution. This situation was explored already in [4] in 
the intermittent regime, but it is also possible to generate homoclons in the 
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Fig. 4. Incoherent Homoclons for ci = C3 = 0.5. The first cohinin shows the evolution 
of \ A\ (white: |^| = 1) while the second column shows the evolution of q (grey: q = 0, 
light: q > 0). The third column shows the evolution of the phase of an artificially 
created "fast" field: for that A has been multiplied with a plane wave solution. (a,b,c) 
Initial phase-twist just below Homoclon equilibrium. (d,e,f) Initial phase-twist just 
above Homoclon equilibrium. (g,h,j) Initial phase- twist further above Homoclon 
equilibrium; the only difference to d,e,f is the shorter lifetime of the incoherent 
Homoclon here due to a different initial condition. 

"laminar" CGLE regime dominated by stable plane waves (see Fig. 4). In ad- 
dition, the profiles of incoherent homoclons evolve, within good approximation 
along a one-dimensionally family; if one takes many snapshots of incoherent 
homoclons and orders them by their extremal phase-gradient qex, profiles with 
the same q^x should look fairly similar (see [4] for a theoretical check of this) . 
These properties open up the possibility for a number of experiments that 
follow below. 

3.2.1 Creating transient Homoclons 

As discussed above, transient Homoclons can be generated in a large range of 
coefficient space, provided one can locally perturb the system sufficiently. A 
possible experimental protocol would be as follows. 

• Generate a stable plane wave, and find a way to locally perturb this wave. As 
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discussed above, thermal systems could be perturbed by local heating and 
many systems could be perturbed by mechanical means. Whatever method 
one uses, one would expect to be able to tune the strength of the perturba- 
tion. 

• Clearly, for small enough perturbations no defects will be formed, but in- 
stead the perturbation will drift with the group velocity and decay diffu- 
sively. The question is: can one make a large enough perturbation that grows 
out to a defect? This is the necessary ingredient; and this is where clever 
experiments will be the best way to find out what perturbations are most 
successful. Prom a theoretical point of view, it is known that perturbations 
that "twist" the phase-field are most effective. I would imagine that the 
effect of local heating is a local increase of the effective value of and the 
resulting change in time and spatial scales will certainly lead to the gener- 
ation of phase-twists. If these are done in the middle of the system, likely 
two pairs of twists are generated (at either side of the perturbed region), so 
it may be advantageous to perform the perturbations at the edge. 

• Denote the control parameter that controls the strength of the perturbation 
as V. If, defects are formed for V = V2, while for V = V\ the perturbation 
decays, the experimental protocol is then to repeat this experiment for a 
number of values between V\ and V2 to obtain a critical value between 
decay and defect generation. Approaching this value, the theory predicts 
that transient Homoclons of longer and longer lifetimes occur, and that their 
lifetime diverges as —\~Hn{}y — V^|) [5], where A is the unstable eigenvalue 
of the Homoclon [5] . 

• The theory predicts that the tuning of a single parameter is sufficient to 
come arbitrarily close to a coherent Homoclon, i.e., if one has two control 
parameters, say V and W ^ than there should be a continuous curve in V, W 
space where the Homoclon lifetime diverges. 

• Finally, if one is able to tune the wavenumber of the plane wave that the 
Homoclon propagates into, this also can be used as a parameter. This prop- 
erty was explored in [5] in the regime where the coefficients C\ and C3 are 
such that the defects generated by Homoclons lead to the birth of new ho- 
moclons and periodic hole-defect states occur. Even when this is not the 
case, i.e., in the "laminar" regime, the lifetime of a Homoclon can still be 
controlled by the asymptotic wavenumber. 

3.2.2 Spontaneous incoherent Homoclon generation 

There are a number of states in which incoherent Homoclons are generated. 
First there is the hole-defect chaos that occurs in a part of the so-called spatial- 
temporal intermittent regime of the CGLE [4,5,19]; as far as I am aware, this 
state has not yet been experimentally observed. In experiments it is often 
possible to have sources between left and right traveling waves, and in that 
case a sufficient lowering of e leads to the creation of unstable sources that 
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send out irregular phase-twists that can lead to transient Homoclons [29] . 

In both of these states one may study the a and q profiles of snapshots of the 
holes, and see if they indeed can be ordered as a one parameter family as is 
the case for the CGLE [4]. In addition, when following q^,^ as a function of 
time for a hole, the time derivative q^x should be a function of qex only [5]; 
this should be revealed in scatter-plots of q^x versus qex- Finally it should be 
remarked that the late stages of defect formation in hole-defect dynamics also 
is expected to show some universal features; for more see [5]. 

3.3 Nozaki-Bekki holes 

I will not spend much time on the incoherent dynamics of Nozaki-Bekki holes 
since there is already a very large body of literature available [7,10,11]. There 
is a potentially large range of dynamical behaviors possible for NB holes, 
due to either dynamical or structural instabilities. When one would observe 
dynamical NB holes in an experiment, it may be interesting to measure the 
time-scales as a function of e. For a CGLE-like structure, these should scale as 
e~^, but if the dynamics is due to the structural instability one may expect a 
larger exponent here, since the deviation from the CGLE can also be expected 
to scale with e; on the other hand, when NB holes are dynamically unstable, 
the structural instability may be irrelevant. 

3.4 Identification of local structures 

What sort of states can be expected to be seen in experiments? In the light 
of the experimental data available and our discussion of the properties of 
the various coherent structures, three possibilities seem most common: (i): 
Stable MAWs (either in small systems or for high winding- number) [21,25,32] 
(a): Transient MAWs that evolve to defects [25,26] (Hi) Holes that may be 
Homoclons or NB Holes. Since cases (i) and (ii) have been discussed already 
in the sections above, only the distinction between Homoclons and NB holes 
will be addressed. 

When one observes coherent holes, this situation appears simple: When the 
two wavenumbers of the adjacent waves are equal, they may be Homoclons, 
when the are different, they may be NB holes. The problem is that in real- 
istic situations, both wavenumbers always will be different, due to noise etc; 
as can be expected due to the structural and dynamical instabilities of these 
structures, typically incoherent holes are observed. For incoherent structures, 
one cannot take a single snapshot and compare this to the profile of either 
Homoclons or NB holes. Also, the wavenumbers at both sides of an incoherent 
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Homoclon are, due to the combination of an evolving core and phase conser- 
vation, no^ equal (see [4,5]). 

There is, however, a clear difference between incoherent Homoclons and NB 
holes, which manifests itself in the direction of the full nonlinear group- velocity 
of the adjacent waves. I would even like to go so far as to propose this as a 
definition: 

In the frame moving with the incoherent hole, incoherent NB holes send out 
waves, while incoherent Homoclons have one incoming and one outgoing wave. 
Needless to say, checks on the direction of the group-velocity are also useful 
for coherent holes. 

An important consequence of this is that two incoherent NB holes cannot sit 
immediately next to each other, but need to be separated by a so-called sink 
[24] (unless they have very different propagation velocities) . This feature can 
be used when it is difficult to measure the group velocity, or when this velocity 
becomes close to the propagation velocities of the holes (which is the case for 
small s, where the group velocity becomes close to the linear group velocity 
Sd, see appendix A). 

I will now briefly discuss some earlier experimental work in which holes are 
observed; for a more general overview of earlier experimental work, the reader 
is referred to [11,15]. In [12], convection of a low Prandtl number fluid in 
an annular container is studied. Heated wires embedded in the walls of this 
container yield horizontal temperature gradients which result in a stable two- 
concentric-roll pattern. These roles become unstable to an oscillatory insta- 
bility when an additional vertical temperature gradient is imposed, leading to 
one-dimensional waves that travel along the anmihis. For Raylcigh numbers 
larger than 1.2 Racru, these waves are stable, and the authors focus on the 
case that Ra — 1.1 Rucru, when the homogeneous wave becomes unstable and 
strong but slowly varying amplitude and wavenumber modulations develop 
along the pattern. The authors conclude that these modulations presumably 
are not described by (coupled) CGLEs; the modulation does, however, yield 
a local area of relative large wavenumber in which holes are generated. 

These holes propagate for a finite amount of time, and in many cases lead to 
the formation of defects. To quote the authors: "We observe the spontaneous 
creation of traveling dips in a traveling wave pattern. They appear in a region 
of low amplitude and propagate through it. They change their depth and the 
value of the phase jump at their core as they mov^' . While it is difficult to 
make definitive statements, it may very well be that the holes are formed 
in a manner similar to that shown in Fig. 3, i.e., due to the instability of a 
wave; the authors in fact claim this in their conclusion. At the point of writing 
of this paper, only NB holes were known and the authors speculate that the 
structures they observe are NB holes, but I doubt that: now that we know that 
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beside NB holes, there are also Homoclons and "ghost" states such as shown 
in Fig. 3, the question what sort of holes where observed in this experiment 
is, in my opinion, open. 

In [13], oscillatory Rayleigh-Benard convection in a low Prandtl number fluid 
(argon gas) is studied. The state of 'coupled oscillators' that the authors ob- 
serve may become chaotic, in which cases holes occur (see their Fig. 6). The 
holes are "characterized by a strong amplitude dip and a phase jump local- 
ized at their cores." They have a finite lifetime and terminate in defects. "The 
slower they move, the deeper the amplitude dip". The incoherent motion of 
the holes, their evolution to defects and the fact that a defect can generate a 
pair of holes with opposite phasejump (see Fig. 6 of [13]) all are strongly rem- 
iniscent of the dynamics of incoherent Homoclons. As before, when this paper 
was written, only NB holes were known and so naturally the authors speculate 
that their holes are of this form, but it seems fair to say that the situation 
is unclear; in fact, the dynamics shown in their Fig. 6 could at least as well 
be due to homoclinic type holes, also because strong evidence for significantly 
different wavenumbers at both sides of the holes is absent [33] . 

More recently, a careful study of holes occurring in hydrothermal waves in a 
laterally heated fiuid layer was carried out. The system consisted of a rectan- 
gular container of length 25cm and width 2cm, in which left and right traveling 
of typical wavelegnth 5mm occured. The coefficients f^^, and tq where all 
measured, and e was at a value of 0.19. The dynamical state that occurs then 
(see Fig. 2 of [14]) is that of an unstable source which separates left and right 
traveling waves; this source sends out perturbations that develop into (incoher- 
ent) holes. The authors claim that these holes are of NB type. By measuring 
the profile of such hole (in particular their left and right wavcnumber and 
the value of the amplitude minimum), and comparing these to values know 
from the NB family, they estimate that Ci —1.5 and C3 0.4; the authors 
claim that similar values of Ci and C3 were obtained for a number of different 
holes send out by the unstable source. Unfortunately, no independent measure- 
ments of the coefficients Ci and C3 could be extracted for the hydrothermal 
waves. Using the values of the coefficients quoted by the authors, 1 studied 
the behavior of the (coupled) CGLEs here, using the values of the coefficients 
quoted by the authors (including a rescaled group- velocity of 1.7), and found 
indeed an unstable source that sends out hole-like perturbations (see [16] for 
an explanation of the instability of the source). In the simulations these holes 
are quite incoherent, and it is difficult to see whether they are homoclinic or 
heteroclinic. 

I have no acces to the raw data, and the evidence for NB holes presented is 
impressive. Nevertheless, there are some puzzling aspects to this experiment. 
For the values of coefficients and wavenumbers quoted, NB-holes are quite 
unstable; I found in simulations of a single CGLE that even the generation of 
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transient unstable NB holes is difficult to achieve for generic initial conditions; 
I do not know of any mechanism that would generate unstable NB holes, and 
it is therefore surprising, but not impossible, that these holes are seen in 
experiment! Another thing that worries me is the size of the holes, and their 
ill separation. The snapshot of a NB hole that the authors show in their Fig. 4 
has an extension of roughly 50mm on the scale of their spacetime diagram 
Fig. 2Q, and so I would have expected to see clear differences in wavenumbers 
and sinks in the spacetime-diagram. But this diagram shows the occurrence 
of holes that are close together (distance approximately 10-20 mm), without 
any indication of a sink in between that would separate them if they would be 
NB holes. Also, NB-holes are sources, so they should send out waves; again, I 
don't see any evidence for this. 

An alternative, but possibly less satisfying interpretation of the data goes as 
follows. First note that the total time shown in the spacetime plot is quite 
short also (105 sec, while ro/e = 8 sec [14]), and that all propagation is 
completely dominated by the linear group velocity. The holes, when they are 
just send out from the source, form defects, and the pre-defect holes appear 
reminiscent of (transient) Homoclons; their dynamics is quite fast and so the 
wavenumbers to the left and right should be substantially different; in fact 
they seem so incoherent, that detailed comparison to either NB or homoclinic 
holes seems fairly hopeless. As far as I understand, the holes studied by the 
authors are the ones formed after these defects have occurred. If one now takes 
the single CGLE and studies defect formation for the coefficients quoted by the 
authors, one finds (for example by taking an initial wave of high wavenumber 
(0.6)), that after the last defect was formed, a transient state persisted for 
a few timeunits that is characterzized by a dip of \A\, and two wavenumber 
"blobs" (that come from the decaying defect); one could interpret these as 
patches of waves with different wavenumbers to the left and right of this dip. 
When one switches on the group velocity, such structures are swept away, 
and could look similar to traveling holes. These structures decays than fairly 
rapidly (of the other of 10-20 time units in CGLE, which would be 80-160sec 
in the experiment), but are, in my opinion, not related to any coherent or 
incoherent hole state. Note that the occurrence of short timescales is inherent 
in small-e experiments in finite systems; either very large systems or periodic 
boundary conditions may make more detailed comparison of the holes send 
out by unstable sources in hydrothermal waves possible in the future, and I'm 
looking forward to see what would come out of such experiments! 

Finally, in [34] , the dynamics of holes and the formation of defects is studied in 
the Taylor-Dean system. In this system a traveling wave with negative phase 
diffusion occurs, and phase gradients concentrate and lead to defects (in re- 
lated work, these authors observed stable MAWs [32]). The dynamics depicted 



^ since ^o/v^ = 2.5 mm [14]. 
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in Fig. 10 and 11 of [34] suggest that transient Homoclon-like structures play 
a role here. Finally, in a recent study of sources and holes in a heated wire 
convection experiment [29] , the dynamics of holes send out by unstable sources 
is studied. The temporal evolution of the profile of the holes shows properties 
reminiscent of that of Homoclons; nevertheless some questions remain open. 

Notwithstanding all these indications, in my opinion neither NB holes or Ho- 
moclons have been observed unambiguously yet. Experiments in which the 
coefficients of the amplitude equation are measured independently and com- 
pared to the profile and evolution of holes, or precise comparison of group- 
velocities and propagation velocities of holes may well show (unstable) NB 
holes. Similarly, experiments in which initial phasetwists can be shown to lead 
to either long or short lived holes (similar to Fig. 4) may yield more definite 
experimental evidence for Homoclons. 



4 Conclusion and Outlook 

In this paper I hope to have given a broad and easy accessible introduction 
to the dynamics of MAWs and Homoclons, and some clarification on the Ho- 
moclon/NB hole issue. Hopefully some of the suggestions made above will 
inspire new experiments. Maybe this is the right moment to give a personal 
perspective on why one could still be interested in complex dynamics in one- 
dimensional systems. 

There is no overall applicable principle or method to describe non-equilibrium 
systems, but a large subclass is characterized by a combination of nonlinearity 
and spatial extend. While in general it is not clear how to go from the tools 
developed for low dimensional dynamical systems to an effective description 
for systems with many degrees of freedom, the coherent/incoherent structures 
framework sketched in this paper seems to be able to form such bridge in the 
case of the CGLE. 

Possibly the greatest advantage of studying one-dimensional systems is that 
their time evolution can be captured in two-dimensional space-time plots, 
which allow the development of intuition for these systems. Without these, 
the discovery of most of the Homoclon and MAW dynamical properties would 
have been much more difficult. 

The uncovered mechanisms by which an unstable wave can give rise to MAWs, 

or, beyond the saddle-node, to defects, and the existence of Homoclons that act 
as separatrices between defect free and defect developing states, are not trivial. 
1 have some hope that these mechanisms may prove to be more general, and 
therefore it will be extremely interesting to see what happens in experiments. 
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Finally, when e is increased sufficiently, new mechanisms will start to play a 
role in the experiments, and the interplay between these and the "low e" dy- 
namics described here will be interesting. The range of experiments, although 
potentially large, over which the CGLE dynamics is applicable has not yet 
been sufficiently mapped out: let us hope that more will be known in the next 
few years. 
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A Coefficients of the CGLE in the labframe 

The spacetime diagrams one obtains in experimental situations cannot be di- 
rectly compared with the predictions from the CGLE, since the simple form 
of the CGLE (1) is only obtained after a number of rescalings and coordinate 
transforms that I will discuss here. It is convenient to introduce first the di- 
mensionless parameter e that measures the distance of the controlparameter 
to threshold [15]. At threshold (e = 0), the mode with dimensional wavenumber 
qc and frequency ujc (nonzero) becomes unstable. 

The central observation underlying the amplitude equation approach is that 
for small but finite positive values of e one expects a hand of wavenumbers 
of width oc ^/e to play a role. Hence close to threshold one expects that a 
physical field u can be written as the product of a slowly varying amplitude A 
with the critical mode: u{xd,td) oc exp{i{qcXd — LJctd))A{x,t) +c.c. [15], where 
X and t denote the slow non-dimensionalized space and time coordinates. Af- 
ter substituting this ansatz into the underlying physical equations of motion 
(assuming that they are known) and performing a rather tedious expansion 
up to third order in s [15,35,36], one then obtains the appropriate amplitude 
equations: the cubic complex Ginzburg Landau equation. In its full dimen- 
sional form, which is most appropriate for highlighting the problems one may 
encounter when comparing to real data, the equation reads: 

dA OA d^A 
To(^ - sa^) = £(1 + ico)Aa + eo(l + ici)^ - go{l - ics)\Aa\''Ad .(A.l) 



23 



To compare experimental data to the CGLE for Ad (Eq. (A.l)), one has to 
ehminate the fast scales corresponding to the critical mode, and to do so one 
proceeds as follows. 

• Perform a Laplace transform on the experimental spacetime data set u to 
obtain a complex valued field C/ as a function of the labframe coordinates 

Xd and td [29]. 

• Find the onset for pattern formation, and measure the critical wavenumber 
Qc and frequency cUc, i.e., characterize the wave obtained for e as close to 
zero as possible. 

• Demodulate the field U to obtain the dimensional field Ad, by writing Ad — 
U exp{i{qcXd — ojctd)) (assuming one wave with positive phase- velocity) . 

One now has obtained a complex field A^ that is described by the dimensional 
amplitude equation Eq. (A.l). 

For a theoretical analysis, most of the coefficients occurring in Eq. (A.l) can 
be scaled out. The coefficients tq and give typical temporal and spatial 
scales to the equation as given by the dispersion relation for plane waves. 
Going to the so-called dimensionless slow scales x and t, which are defined 
as X := ^/e^Q^Xd and t := ET^^td sets the coefficients {o,t"o and the linear 

— 1/2 

growthrate e equal to 1. By writing Ad = Qq exp (icot)A, the coefficients qq 
and Co are scaled to 1 and respectively, and the non-dimensionalized form 
of A is obtained. The linear, dimensional group velocity Sd can be removed 
by going to a comoving frame, but this only makes sense when the system 
has periodic boundary conditions; for a finite system with fixed boundaries, 
the value of Sd does play an important role. Assuming that one can go to the 
comoving frame, the CGLE in its simple form (1) is obtained. 



A.l Measuring the CGLE coefficients 

The coefficients tq, Sd, cq, .^g, Ci, (jfo and C3 can in principle be calculated from 
the underlying equations of motion [15,35,36]. For many experimental situ- 
ations, however, such calculations are not available or can only be done in 
certain approximations; for some systems the equations of motion or bound- 
ary conditions arc not even known in enough detail to allow such calculations. 
Since the scale and even qualitative character of the dynamics depends on 
these coefficients, their knowledge is essential. The task of measuring these 
coefficients appears rather nightmarish at first sight, but as I hope to show 
below, may in fact not be terribly complicated. For recent examples of exper- 
imental determinations of these coefficients see [14,29,26]. 
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A. 1.1 Homogeneous solutions and dispersion relation 

When one searches for homogeneous but possibly growing or decaying plane 
wave solutions of Eq. (A.l) of the form a(td) exp{i{qdXx — uj(itd))\^oiae obtains 
a complex valued equation, of which the real part reads: 

da/dtd = e- ^Iql - g^a^ . (A. 2) 



When a is time- independent, the real and imaginary part of this equation are: 



Touj = eco - TQSdQd + ilciq^ - goc^a^ (A. 3) 

= ^-eod-^?oa' • (A.4) 

Finally, when one restricts oneself to plane wave solutions, where Eq. (A.4) 
is satisfied, Eq. (A. 3) becomes the "nonlinear" dispersion relation for plane 
waves: 

tquj = e{co - C3) - ToSdqd + ^o(ci + cs)?^ (A. 5) 



These equations will be the basis for the experiments described below. 



A. 1.2 Quenches of e 

The coefficients tq, cq and C3 can be obtained by performing experiments where 
e is suddenly changed from one value to another. The simplest case is when 
initially the system is below threshold {e < 0) and then suddenly the control 
parameter is changed to a positive value A nonlinear plane wave state will 
start to grow then. Assuming that one is in the linear regime, that this growing 
wave is spatially homogeneous and has wavenumber = 0, one finds from 
Eqs. (A. 2) that a oc exp{TQetd). Denote the time interval during which the 
wave goes from 10% to 50% of its final strength by ti. Plotting e versus 1/ti 
then should show a linear relationship with a slope give by tq. 

During such quenches one can measure the frequency u, both when the wave 
just starts to grow and when it is fully developed. From Eq. (A. 3) one finds 
that, for Qd = 0, an infinitesimal wave has frequency e/tqCq, while a fully 
developed wave has frequency (Eq. (A. 5) £/ro(co — C3). So from these two 
measurements Cq and C3 can be determined. 

Of course, the assumption that the waves are homogeneous and have wavenum- 
ber Qd = may not always be satisfied. However, it is known that for small 

^ Note that the critical wavenumber and frequency have already been split off: the 
bare labframe wavenumber qi^b is equal to Qc + Qd- 
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but positive e the band of allowed wavenumbers shrinks oc while the spa- 
tial modulational scales similarly. Therefore, when the growth-rate is rapidly 
changed from Si to £2 and back a number of times, where Si is close to zero 
while £2 is not, one expects to indeed have a homogeneous plane wave of 
wavenumber very close to zero. Then one can use Eq. (A. 2) which, for q^i — 0, 
becomes da/dtd — e — goo^ to get a full numerical prediction of a as a function 
of time. Comparing this full solution to the experimentally obtained curves 
for a{td)-i both for the "up" as well as for the "down" quench for a range of 
values of £1 and £2 then gives accurate estimates for both qq and tq. 

A. 1.3 Propagation of linear perturbations 

Suppose one has generated a stable plane wave and locally perturbs this wave. 
Then, as discussed above, the temporal evolution of this perturbation will be 
a combination of slow diffusion and advcction with group velocity dujd/dqd- 
For waves of wavenumber = 0, this group velocity is the linear group veloc- 
ity Sa, and for more general waves one finds by differentiating the nonlinear 
dispersion relation that the nonlinear groupvelocity is sa + 2^Qq(i{ci + c^/tq. 
For a measurement of Sd it is therefore easiest to make qd equal to zero, which 
can, for example, be done by varying e up and down as described above. It 
should be noted that while small perturbations may be difficult to observe in 
noisy snapshots of the system, they often become quite clear in space-time 
plots of the raw data u. 

When the perturbations are relatively smooth, such that da/dt is small, a 
comparison of the instantaneous amplitude a and local wavenumber q leads to 
another condition on the coefficients of the CGLE. From Eq. (A. 3) one finds 
grQ|>l|2 — £ _ ^^q'^. Plotting |Ap/£ versus q'^/e one expects a hnear relation. 
For I the value of |>lp/£ approaches go, while the slope of the curve will 
be equal to Co/9o- 

A. 1.4 Fronts 

When e is suddenly increased from below to above threshold, details of the 
noise and the initial conditions determine the details of the growth of the 
nonlinear state. Above it has been assumed that one has a homogeneous state, 
but when a localized perturbation is applied just before such a change of e 
occurs, it is possible to create a pair of fronts that propagate into the unstable 
a = state, and leave a nonlinear state in their wake. The point is that for 
large times the velocity of the front and the wavenumber of the nonlinear state 
are selected and can be calculated from an essentially hnear analysis [15,24]. 



These velocities are Sd ± 2^q/tqJ e{l -\- cf), and the wavenumber of the plane 



wave they leave behind is qc -\- J e / {1 -\- c\) [24]. It should be noted that 
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the velocity and waveumber only relax to their asymptotic values as so it 
is hard to avoid making some errors here. 



A. 1.5 Eckhaus instability 

As discussed in the section on MAWs, one can manipulate the effective wavenum- 
ber of a plane wave state by changes of e. At some point, when the effective 
value of q becomes too large, one will encounter the Eckhaus instability, which 
for general e is given by < ^CcT^sZ^^^f^ [20]. This could be used as an 
independent check on the coefficients. 



B Coherent structure ODEs 



The coherent structure ODEs obtained when the ansatz Eq. (4) is substituted 
into the CGLE can be written in a number of forms and for a number of 
different dependent variables. All useful representations will have a and q :— 
9^0 as variables, but there are broadly speaking two choices; either b :— d^a or 
K := {d^a)/a [24,4,2,3]. The latter representation is particularly useful if one 
wants to study structures, such as fronts, that asymptotically decay to A = 0. 
K then measures the exponential decay of the profile to the zero state. Using 
this latter notation the set of coupled ordinary differential equations (ODE's) 
becomes 



d^a = na , 



1 



1 + ici 



-1 — iuj + {1 — ic3)a^ 



vz 



(B.l) 
(B.2) 



where the complex quantity z is equal to K + iq, which is also equal to \n{A). 
Equation (B.2) is equivalent to two real valued equations, so (B.l, B.2) can also 
be seen as a 3D real- valued dynamical system [24] . Note that since the CGLE 
is a complex, second order equation one may have expected to find 4 coupled 
ODEs, but they can be reduced to 3 equations using the phase symmetry of 
the CGLE. 



For completeness note that another ansatz is frequently encountered in the 
literature: A — a(^) exp{i{qx —ujt)) exp(i0(^)). This may be convenient if one 
wants to split off an explicit wavenumber q and require that d^(f) goes to zero 
for ^ — ±oo. It is easy to show that this ansatz is equivalent to ansatz (4), 
and so the apparent three free parameters of this ansatz can be reduced to 
two [11]. 
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